Analytic estimates for penguin operators 
in quenched QCD 
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Abstract 

» . 

Strong penguin operators are singlets under the right-handed flavor symme- 



try group SU(3)^. However, they do not remain singlets when the operator is 
embedded in (partially) quenched QCD, but instead they become linear com- 
binations of two operators with different transformation properties under the 
(partially) quenched symmetry group. This is an artifact of the quenched ap- 
proximation. Each of these two operators is represented by a different set of 
$— i ' low-energy constants in the chiral effective theory. In this paper, we give ana- 

lytic estimates for the leading low-energy constants, in quenched and partially 
quenched QCD. We conclude that the effects of quenching on Qq are large. 
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Even at present, the number of light dynamical (or sea) quarks in Lattice QCD 
computations is often not equal to that of the real world. In particular, in quite a few 
of the more difficult computations, no sea quarks have been taken into account at all 
- in other words, they have been done in the quenched approximation. 

The adaptation of chiral perturbation theory (ChPT) to the quenched theory has 
been useful for gaining insight into the effects of quenching. However, ChPT itself 
does not give any insight into the values of the parameters of the effective theory, 
the so-called low-energy constants (LECs). It follows that ChPT also does not tell us 
anything about the effects of quenching on the values of the LECs. Within Lattice 
QCD, therefore, we do in principle not know anything systematic about the errors due 
to quenching, until computations are done with dynamical quarks, and the results are 
compared with those of quenched QCD. 

At the same time, analytic approaches are being explored to obtain estimates for, 
in particular, electroweak-interaction LECs. Some of these approaches make sophisti- 
cated use of analytic knowledge available about QCD, such as its chiral behavior, the 
operator-product expansion, and large- N c techniques. For the purpose of this paper, a 
relevant reference is Ref. [T], to which we also refer for references to other work. While 
some assumptions have to be made in such approaches because a non-perturbative 
analytic solution to QCD is not available, they are often so tightly constrained that 
it is quite likely that the results obtained from them will help us with our quantita- 
tive understanding of hadron phenomenology. Moreover, they have the advantage of 
exhibiting the underlying reasons for the size of specific contributions. 

It is therefore natural to adapt these analytic techniques to include the effects of 
(partial) quenching. This is useful, because it gives us quantitative information about 
the effects of quenching. Results from this approach can thus serve as a guide to the 
merits and pitfalls of the quenched approximation. 

In this note, we apply the large- iV c expansion to the calculation of the leading LECs 
associated with the strong penguin Q 6 [2] - It was recently pointed out that even the 
definition of the quenched version of this operator is ambiguous and that as a 
consequence, the quenched theory has more LECs associated with this operator than 
the unquenched theory. It turns out that the leading-order LECs can be estimated 
analytically if one ignores order 1/N% corrections. These estimates tell us that the 
effects of quenching on the contribution of Qq to non-leptonic kaon decays are likely 
to be large. 

The operator Q e is defined as 

Q 6 = 4(^7X) E (3S7m«S). (!) 

q=u,d,s 

where q^L = Pr,lQ with P R L = |(1±7 5 ), and a and (3 are color indices. This operator 
transforms in the (8, 1) representation of SU(3)lxSU(3)#. 

In order to "embed" this operator in the quenched theory, let us briefly recall how 
one may define quenched QCD as a field theory jH|3]. For each quark q one introduces a 
ghost quark q with the same mass, spin, and color, but opposite statistics. The opposite 
statistics cause the path integral over the ghost quarks to cancel the quark determinant 
(for each gauge field configuration), effectively replacing the quark determinant by one. 



1 



This is precisely the definition of the quenched approximation. It follows that the flavor 
symmetry of the quenched theory is not described by SU(3), but by the larger, graded 
group SU(3|3). Qq is not a singlet under SU(3|3)#, but instead can be decomposed as 


q 6 = \qT + q2 ns , (2) 

i> 

Q Q , NS = 4(32tX) £@*7^S) , 

N = ^diag (1,1, 1,-1, -1,-1), 

where ip = (q,q). Qq S clearly is a singlet under SU(3|3)#, but Q$ NS is not. 

To leading order in quenched ChPT, these operators are bosonized by (we work in 
euclidean space) 

Qf -> -ag'^str (AL M L M ) + ag'^str (25 A(SM + MS f )) , (3) 
q«*s _> / 2 «f s str (ASiVSt) . 

In these expressions £ = exp(2?$//) is the non-linear field describing the quenched 
Goldstone-meson multiplet, M is the quark mass matrix, Bq is defined in Ref. pi . and 
/ is the chiral limit of the pion decay constant normalized such that f n = 132 MeV. 
A is the tensor picking out the octet operator sd = qAq, and L M = iHd^T) is the left- 
handed current to leading order in ChPT. the LECs associated with the order 

p 2 weak octet kinetic and mass terms [7., while a^ s is the order p° LEC associated 
with the non-singlet operator Qq Note that, while this operator is enhanced in 

ChPT relative to the singlet operator, its contributions to matrix elements with only 
physical quarks on the external lines start only at order p 2 IHj- The subscript q is 
there as a reminder that the a's refer to LECs of the quenched theory. 

We begin with estimating the magnitude of a^ s . A simple way of doing this was 
outlined in Ref. One first rotates Qq NS by an SU(3|3)^ rotation into 

Q& NS = 4(S£V?) E(^7m^S) , (4) 

= -4 {(sP R q) (qP L d) + (sP R q) (|P L d)) 

where in the second line we fierzed the operator, paying careful attention to the fact 
that the ghost-quark fields are commuting. One then considers the K° — > matrix 
element of this operator, with K° a hybrid kaon made of a ghost-o? quark and a physical 
anti-s quark. The advantage of considering this matrix element is that it is of order p° 
in ChPT 0: 

(0\Qr S \K°) = 2 i fa» s + O(p 2 ). (5) 
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It turns out to be quite simple to find an expression for a^ s . By carrying out Wick 
contractions, we may write the K° — > matrix element as 

(0|g 6 QiVS |^ ) = ((ss}-(dd})(0\s l5 d\K }(l + o(^)) (6) 



-4 (0 1 (sP R qq P L d + sP R fq P L d) \ K°) 



where qq denotes the contraction of q with q, and likewise for qq. The first line contains 
terms with two quark loops, while the second term has only one quark loop. In order 
to connect the two loops on the first line, at least two gluon lines are needed to obtain 
a non-zero contribution, making such connected contributions suppressed by 1/iVp. To 
leading order in 1/N C we thus obtain the factorized contributions shown explicitly. 

The one-loop contribution maybe dealt with as follows. First, we observe that ghost 
and physical quark propagators are equal, flavor by flavor, by construction. Then, we 
may, in this term, rotate the K° back to a K°, and correspondingly the di quark to 
di. For the second term we find 

-A(0\(sP R qqP L d + sP R ^_P L d)\k } = -8(0\sP R qqP L d\K°) (7) 

= -8(0\(sP R q)(qP L d)\K°) . 

The last step follows because, as one can show, the Wick contractions leading to contri- 
butions with two quark loops cancel each other in the chiral limit. The last expression 
in Eq. (JZj) is just the K° — > matrix element of Q§ 2 which is of order p 2 . This is also 
true in the quenched theory [Hj. Hence, this term does not contribute to af? . 

Using that, in the chiral limit, (dd) = —(ss) = ^f 2 B , we thus find for the matrix 
element to order p° that 

(0\Q% NS \K°) = 2(ss)(0\s l5 d\K°) (8) 

= 2{ss)if^^ = -ifBl 
m s + m d 

Comparing with Eq. (jSj), we obtain an estimate for a^ s correct to order 1/N 2 : 

< S = -\f 2 Bl (9) 
In order to get an idea about the value of we may compare it to the value 

(8 1) 

of a i . Fortunately, it turns out to be remarkably simple to obtain an estimate for 

a q { in the quenched theory. It turns out that in the quenched case, the unfactorized 
contribution vanishes! This can be seen as follows. The fierzed form of the singlet 
operator Q^ S is, from Eq. (J2J), 

Qf = -8 ((sP R q)(qP L d) - (sP R q)$P L d)) . (10) 



2 through the weak mass term 0. 
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Following the analysis of Ref. pQ, 3 again there is a contribution with two quark loops, 
and a contribution with one quark loop. The one-loop contribution again corresponds 
to the terms in which and q R or q R an d q R are contracted, and thus vanishes because 
of the relative minus sign in Eq. (JTUJl . This leaves us with the two-loop contribution, 
which factorizes to order 1/N 2 , yielding [5| 

ag' 1} = -8L 5 f 2 B 2 (11) 

(equivalent to g$ = —1QL 5 B 2 /F^, F = f/y/2 in the notation of Ref. [JJ). 4 Note that 
quenched L 5 does not run. 

Putting things together, we find that 



a NS 



a ql 



(8,1) 




(12) 



Our first conclusion based on these results is that a™ is likely to be large compared 
to the singlet LECs. The value of L5 is of order 10 -3 both in the quenched [HI] and 
unquenched |2J theories, making this ratio of order 60. This is not small, 5 and casts 
doubt on the tentative conclusion of Ref. [TT] (concluding section) on the size of a^ s . 

If otq S is not small, this could have a dramatic effect on the extraction of from 

K° — ► 0, which in turn is needed for the extraction of a^' 1 ^ from K — > rr [T_, because 

it appears at the same order as o& in ChPT in the K° — > matrix element [3]. 

' (8 ~L) 

Another interesting observation is that the value of a q { maybe significantly smaller 
(in absolute value) than that of the unquenched theory. This is because of the absence 
of unfactorized contributions in the quenched theory, which, in the unquenched theory, 
have been estimated to be sizeable compared to the factorized contribution, and of the 
same sign [JJ. 6 Even if one would decide that one should define Qq in the quenched 
theory to be Q® s only, omitting altogether, as an alternative possibility [31 El/ 

(8 1) 

this smaller value of oc q { would lead to a reduction of the value of e'/e in quenched 
QCD, relative to its unquenched value. 

It is rather easy to extend these estimates to the partially quenched situation, in 
which N (massless) sea quarks are added to the quenched theory JTH]. First, the non- 
singlet part of Qq is now represented by the order p° LEC a^ 8 ' 8 ^ [3J, because in the 
partially quenched case the non-singlet operator is in the same irreducible representa- 
tion as Qs [2] (and thus it corresponds to g ew of Ref. [JJ). One finds that the same 
expression as given for a^ s in Eq. (jSJ) is also valid in the partially quenched theory 

for c^ 8,8 ). For at] 8 , a naive estimate would be to interpolate linearly in the number 

3 As in Ref. p^, we will work in the leading- log approximation, in which Qg does not mix with 
any other operator in the quenched theory. 

4 ^5 is one of the Gasser-Leutwyler constants 6 . 

5 In particular, "strategy 3" of Ref. [S| would yield results significantly different from the other 
strategies. 

6 The unfactorized contribution is found to be about twice the factorized one, at the scale of the p 
mass. 

7 This change is known to have a large effect El EH • 
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of sea-quark flavors N between the quenched and unquenched theories. Using the re- 
sults obtained in Ref. |T], which considers the unfactorized contribution proportional 
to the number of light flavors (which is three in the real world), this would lead to the 
estimate 

a (8,i) _ _g£ 5 y 2 £?2 _| x ( un f ac t or i zec l contribution of Ref. [T]) , (13) 

3 

where of course now L 5 , /, B , etc. take their values in the partially quenched theory 
with N sea quarks. 

Finally, let us comment on the fact that, since our aim was to extract values for 
the leading order LEC for each of the operators we considered, all calculations were 
done in the chiral limit. It is well known that the chiral limit of the quenched theory is 
hampered by severe infrared divergences, and, in fact, the chiral limit of the quenched 
theory may not exist [HlEIlEj- However, we believe that our results for LECs are not 
affected by this issue. If the quenched effective theory makes any sense, its parameters, 
which are the LECs, should be well defined and finite in the chiral limit. All the sick- 
nesses associated with the quenched infrared behavior should be correctly reproduced 
if the appropriate fields, in particular the rj', are kept in the effective theory. 
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